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ABSTRACT 

We examine the usefulness of the strong lensing statistics to constrain the evolution 
of the number density of lensing galaxies by adopting the values of the cosmological 
parameters determined by recent WMAP observation. For this purpose, we employ 
the lens-redshift test proposed by Kochanek (1992) and constrain the parameters 
in two evolution models, simple power-law model characterized by the power law 
indexes v n and v v and the evolution model by Mitchell et al. (2005) based on CDM 
structure formation scenario. We use the well-defined lens sample from the Sloan 
Digital Sky Survey (SDSS) and this is similarly sized samples used in the previous 
studies. Furthermore, we adopt the velocity dispersion function of early-type galaxies 
based on SDSS DR1 and DR5. It turns out that the indexes of power-law model are 
consistent with the previous studies, thus our results indicate the mild evolution in the 
number and velocity dispersion of early- type galaxies out to z = 1. However we found 
that the values for p and q used by Mitchell et al. are inconsistent with the presently 
available observational data. More complete sample is necessary to withdraw more 
realistic determination on these parameters. 
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1 INTRODUCTION 

Since the pioneering work by Turner, Ostriker and 
Gott (Turner, Ostriker & Gott 1984), the strong 
lensing statistics is known as a powerful tool to con- 
strain the cosmological parameters, especially the 
cosmological constant (Fukugita, Futamase, Kasai & 
Turner 1992), as well as to investigate the evolu- 
tion of number density of lensi ng galaxies (iMao 1991 , 

| Mao fe Kochanek 19941 . lOfek. Rix fc Maoz 2003 . 

Mitchell et al. 20051. IChae. Mao fe Kang 2006[ and 
ICapelo &; Nataraian 20071 ). However this method is 
not competitive with other approaches such as Cosmic 
Microwave Backgr ound(CMB). In fa ct recent observation 
by WMAP(e. g. ISoergel et al. 20071 ) determined the cos- 
mological parameters within the accuracy of a few percent. 
The reason is mainly because of the shortage of observed 
lens systems as well as the partial knowledge of the redshift 
information of the sources and lensing galaxies. 

Specifically, the current largest complete lens sample 
comes from the Cosmic Lens All S ky Survey (CLASS; 
iMvers et al. 20031 . iBrowne et al. 2003h , an extension of the 
earlier Jodrell Bank/ Very Large Array Astrometric Sur- 
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vey (JVAS; IPatnaik et al. 19931 . IPantaik et al. 19921 and 
lAugusto et al. 20"oih , but the number of lensing events are 
only 13 lensed radio sources. On the other hand, the SDSS 
Quasar Lens Search (SQLS) has already discovered approx- 
imately 20 new strongly lensed quasars, so this sample is 
expected to become the largest statistical sample of strongly 
lensed quasars. 

On the other hand, there have been no well motivated 
and established theoretical models for the evolution of the 
number density of the lensing galaxies, and there seems 
no definite observational conclusion if the number density 
evolves or not. 

Thus it will be timely to investigates the evolution of the 
number density of the lensing galaxies by adopting the val- 
ues for the cosmological parameters determined by WMAP 
and using a new sample based on SQLS. 

In this paper, we focus on two typical redshift evolution 
models: One model is described by the power-law evolution 
for the number density and the characteristic velocity dis- 
persion as (1 + z) Un and (1 + z) Uv , respectively, (see Chae 
& Mao (2003) in detail). The other is the model based on 
CDM structure formation scenario used by Mitchell et al. 
(2005) which we call Mitchell's model. We use 'lens-redshift 
test' proposed by Kochanek (see Kochanek (1992)) to cal- 
culate the likelihood as a function of the parameters in the 
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models. This method has advantage not to use the total 
lensing probability and thus not to require the knowledge of 
the magnification bias in the sample. We also use the SQLS 
data whose separation is under 2 arcsec as the lens sample 
in order to select only the events caused by a single galaxy 
which is known to be well modeled by Singular Isothermal 
Sphere(SIS) or Singular Isothermal Ellipsoid (SIE). We have 
used the parameters in the velocity dispersion function given 
by SDSS DR1 and DR5. Thus we avoid ambiguities of the 
lens model as possible as we can. 

The organization of this paper is as follows. In section 
2, we explain some basic of lens statistics and show the lens 
models and the evolution models we adopt. In section 3, we 
describe the data that we use. We calculate maximum like- 
lihood value and confidence level of the parameters of each 
evolution model using the differential lensing probability. 
In section 4, we present the results. Finally some discus- 
sion is given in section 5. Throughout this paper we adopt 
A dominant totally flat cosmology with the matter density 
On 0.27. 

2 BASIC OF STRONG LENSING STATISTICS 

2.1 The optical depth 

Following to Turner, Ostriker & Gott (1984), we consider 
the light ray propagating from source (QSO) to observer. 
The differential probability dr of the ray encountering a lens 
in traveling the path of the redshift interval dzL is given by 

cdt 

dr = UL^ZLjVL-, — dzLB (1) 
dzL 

where ul(zl) is number density of the lens (see 2.3), o~l is 
the lensi ng cross section (see!2.2p, and B is the magnification 
bias fsee Turner. Ostriker fe Gott 19841 ). The quantity j^- 
is calculated in the FLRW geometry to be 

cdt _ Rp 1 

dz L ~ i + z L ^/n m (i + ^ L )3 + n A 

where Ro = 1/Ho is Hubble distance, Q m is total mass den- 
sity of the universe, Qa = Ac 2 /3Hq is the normalized cos- 
mological constant, and zl is lens redshift. 

Using JT}, the differential optical depth of lensing in 
traversing dzL with angular separation between AO and 
AO ± dAO is calculated as d 2 r/dz L dA0. 



CTL2im — 7T\/l — e 2 E — O LUm (4) 

where e represents ellipticity and this ellipticity is con- 
strained asO ^ e < 1/5. The SIS model is recovered in 
the e = for aL2im case and this upper constrain e < 1/5 
comes from the fact that the density contours must be con- 
vex, which is reasonable for an isolated relaxed system. And 
Be is Einstein radius for SIS lens and can be written as 

where a is one-dimensional velocity dispersion and Dls and 
Dos are the angular diameter distance between the lens and 
the source, the observer and the source, respectively. 

2.3 Number density of galaxy lens 

The number density of galaxies with velocity dispersion 
(we call this function velocity function after this) lying be- 
tween a and a ± da is assumed to be d escribed by modified 
Schechter function (|Sheth et al. 200l iMitchell et al. 20051 ) 
of the form, 

where 0* is the integrated number density of galaxies, a* is 
the characteristic velocity dispersion, a is the low- velocity 
power-low index and (5 is the high- velocity exponential cut- 
off index. These parameters in the velocity function, (/>*, a, [3 
and a*, are determined by observations. T(x) represents the 
gamma function. The present day comoving number density 
of galaxies can be calculated by integrating (J6j) . 

For these parameters, we use the measured veloc- 
ity dispersion funct ion from SP SS DR1 dSheth et al. 20031 : 
IMitchell et al. 20051 ) and DR5 (|Choi fe Vogelev 20071 ) for 
early- type galaxies. The measured values that Chae (2007) 
outlined are as follows, 

(</>*,a*,a,l3)DRi = [(4.1 ±0.3) x 10" 3 h 3 Mpc -3 , 
88.8 ±17.7 km s~\ 
6.5 ±1.0, 1.93 ±0.22] (7) 

(0*, a*, a, f3) DRb = [8.0 x 10~ 3 h 3 Mpc" 3 , 
161 ± 5 km s~ , 

2.32 ±0.10, 2.67 ±0.07] (8) 
where h is the Hubble constant in units of 100 km/s Mpc -3 . 



2.2 The lens models 

We use a singular isother mal sphere (SIS) a nd a singu- 
lar isothermal ellipsoid (SIE) (Asada et al. 20031 ) as our lens 
model. 

The SIS model is widely used, so we don't explain its p rop- 
erty, (see, for exa mple, iTurner, Ostriker &; Gott 1984 and 
iMaoz fe Rix 19931 ). 

The SIE model can be solved analytically (Asada et al. 
2003). We refer the readers to Asada et al. (2003) for specific 
calculation. We only show the cross sections of this model 
for 4images ((JLUm) and 2images (aL2im) here. 

3tt 4e 2 2 /„x 

CTLAim — — r- =t> E [o) 

o vl — e 



2.4 The lens redshift distribution 

In this subsection we show two evolution models of number 
density of galaxies for J6]). 

2.4-1 The power-law evolution model 

The first evolution model is a power- law type a s most au- 
thors adopted (for example, IChae &; Mao 20031 ). In this 
model, we allow the number density of galaxies 0* and the 
velocity dispersion a* in JB} to vary with redshift as 

=Ml + z)" n (9) 
a*(z) =a*(l + z) Uv (10) 
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The no evolution model corresponds to (y n — 0, v v — 0). 
Note that the case (y n > 0, v v < 0) is consistent with the 
hierarchical model of galaxy formation with bottom-up as- 
sembly of structure. 



2.4-2 Mitchell's evolution model 

The next evolution model is proposed by Mitchell et al. 
(2005). We briefly explain this model. They have taken into 
consideration the number density evolutions by combining a 
halo mass function with a velocity dispersion as a function of 
mass and redshift. They use the prediction by the extended 
Press-Schechter model of structu re formation, which is cal- 
ibrated by N-body simulations ([Sheth &; Tormen 19991 ) as 
the halo mass function. This halo mass function at epoch z 
is introduced by Sheth & Tormen (1999), 



1/2 



exp 



-qv 
2 
(11) 

where p is the mean density, v(z) = 5 2 c / Vf(iW», 6 C = 1.686 
is the extrapolated linear overdensity of a spherical top-hat 
perturbation at the time it collapse, and a 2 (M,z) is the 
variance of the density field at epoch z in linear perturbation 
theory, smoothed with a top hat filter of radius R(M) = 
(3M/47rp) 1//3 . The parameters, (p, q, A(p)), are given by N- 
body simulations in cold dark matter models as p — 0.3, q = 
0.75, A(p) = 0.3222 fsee lSheth fe Tormen 19991 ). 

Next, we sho w the velocity dispersion 
(|Newman fe Davis 200(51 : iBrvan fe Norman 1998h used 
in Mitchell et al. (2005), 



a(M, z) = 92.3A vir (z f ) 1/6 E(z f ) 1/3 ( 1Q13 ^- 1M J km s" 1 

(12) 



M 



1/3 



where 



A vir (z) = 18tt 2 + S2[Q(z) - 1] - 39[Q(z) - l] 2 



and 



Q(z) = 



0+f) 
E*(z) 



(13) 



(14) 



where E{z) = H(z)/H = [H m (l + zf + Q A } 1/2 and z f is 
the mean for mation redshift for a hal o mass M observed at 
redshift z ( see [Mitchell et al. 20051 and lLacev fc Cole 19941 ). 
Combining equation (|11[) and (|12|) yields the velocity func- 
tion at redshift function. 

Following above mentioned procedure, we can compute 
the ratio of the velocity dispersion function at two epochs, 
4>(cr; z)/<f>(a\ 0), as a function of cosmological parameters. 
This ratio can be combined with the velocity function J6}. 



from a survey that is complete and has well- characterized, 
homogeneous selection criteria. From CASTLES data base, 
we choose 56 samples discovered (re-discovered) by SDSS 
from CASTLES data. To eliminate biases in terms of the 
redshift of the lensing galaxy, we have selected systems with 
0.6 < z s < 2.2. This source redshift range comes from SDSS 
selection criterion of the statistical sample of strongly lensed 
quasars (Inada et al. 2007b). 

Furthermore we exclude the systems whose lensing 
galaxies are the spiral galaxies because the lens model of 
spiral galaxies is not well established, and the fraction of 
spiral galaxies in the lensing galaxies has not been decided 
yet. Furthermore the contribution of the spiral galaxies on 
the lensing probability is expected to be small. We also ex- 
clude the systems whose image separation are large (> 2.0 ) 
or the systems that contain multiple lensing galaxies. 

Using above procedure, our sample contains the lensing 
systems whose lens redshift is unknown. We exclude these 
lensing systems because we need the lens redshift and the 
source redshift information in using this lens-redshift test. 
Lenses at a low redshift are more likely to have a spectro- 
scopically confirmed redshift than lenses at high redshift so 
this may bias the sample to systems where the lens is at an 
artificially low redshift. But our aim is to investigate the use- 
fulness of the lens-redshift test in using the lensing sample 
from SDSS only. For this reason, we simply neglect these 
event at this time. In the end our lensing sample consists 
only of the systems with single early-type galaxy as lensing 
object and this contains 13 lensing systems (see Table 1). 

3.2 The maximum likelihood method 

We use a maximum likelihood estimator in our statistical 
analysis of lens redshift distributions. In this paper, we fol- 
low the method proposed by Kochanek (1992) and use the 
differential optical depth to lensing with respect to the an- 
gular separation AO as the probability density. For each sys- 
tems, we calculate the probability density P(zl\zs, AO, X) 
which is normalized unity, where X is the set of galaxy evo- 
lution model parameters. For power-law evolution model, 
this set is X = (z/ n , v v ), and for Mitchell's evolution model, 
X = (p, q). The probability distributions of lens redshifts 
zl with a given angular separation are calculated as follows, 

v>( 1 a, ^ gjfgg^ z ± A9 ' x) 

P{z L \z s , A$, X) = 3— — — — . (15) 

dSS\. z s, AO, X) 

The likelihood as a function of X is given by, 



N 



L{X) = Y[Pi(z L \z s , AO, X), 



(16) 



where N is total system number. 



3 DATA AND METHOD 
3.1 The lens galaxy sample 

Our lens sample is primarily drawn from CASTLES (Muhoz 
et al. 1998) data bas^E While this total sample size of CAS- 
TLES is 100 systems, a statistical analysis requires a sample 

1 http:/ /cfa- www. harvard.edu/castles/index. html 



4 RESULTS 

We investigate two evolution models by estimating the like- 
lihood function (fT6|) . 

First, we obtain constrains on the power- law evolution 
model parameters, (y n , v v ). We calculate the v n -v v contours 
using the SIS lens and the SIE lens for SDSS DR1 and DR5 
data. The results are shown in Fig. 1. The maximum L value 
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Table 1. The columns are: (1) lens number; (2) lens name; (3) source redshift, z s \ (4) lens redshift, z\\ (5) lens separation, A#, in 
arcseconds; (6) grade for the likelihood that the object is a lens: A=I'd bet my life, B=I'd bet your life, and C=Fd bet your life and you 
should worry (CASTLES); (7) Number of images corresponding to each source component, E means extended and R means there is an 
Einstein ring (CASTLES); (8) references. 
List of references: 

1 - Inada et al. (2006); 2 - Lubin et al. (2000); 3 - Ofek et al. (2006); 4 - Eigenbrod et al. (2006a); 5 - Rusin et al. (2003); 6 -Schechter et 
al. (1998); 7 - Fassnacht & Cohen (1998); 8 - Myers et al. (1999); 9 - Eigenbrod et al. (2006b); 10 - Morokuma et al. (2007); 11 - Inada 
et al. (2006); 12 - Inada et al. (2007a); 13 - Burud et al. (2002); 14 - Ratnatunga et al. (1999); 15 - Lehdr et al. (1993); 16 - Treu & 
Koopmans (2003); 



and the la limits are (y n — — 2.0^3 0, v v — O.ISIq - ^) and 
(v n = — 1.0±|-J, v v = 0.0t°-j°) using the SIS lens for SDSS 
DRl(top left panel in Fig. 1) data and DR5 data(top right 
panel in Fig. 1), respectively. In case of SIE lens model, we 
have (y n = l-5±l:g, v v = -0.15lg;?g) and (v n = 2.0±i%, 
Vv = -0.25±S".2o) for SDSS DR1 (bottom left panel in Fig. 
1) data and DR5 data(bottom right panel in Fig. 1), respec- 
tively. 

We compare the estimation of parameters used in the 
power- low evolution model with previous results. Ofek, Rix 
& Maoz (2003) give the la limit dlogcr*/dz = -0.1 
for the evolution of the characteristic velocity dispersion. 
The la limit of the characteristic velocity dispersion of the 
lenses at z = 1 must be within 69%-91% of the present value. 
Meanwhile, Chae & Mao (2003) obtain limit v v = 0.2±°;g 
(la) without using the lensing rate. The la limit of the char- 
acteristic velocity dispersion of the lenses at z = 1 must be 
within 57%- 174% of the present value. On the other hand, 
our results are 84%-104% (the SIE lens for DR1) and 73%- 
97% (the SIE lens for DR5). Thus we obtain the constrains 
on the evolution of the characteristic velocity dispersion 
which are consistent with these two previous studies. 

Ofek, Rix & Maoz (2003) also give the la limit 
dlog 0*/dz = +0.7l} - 2 f° r ^ ne number density. The la limit 
of number density of the lenses at z = 1 must be 30% of the 
present value. Our results are 35% (the SIE lens for DR1) 
and 50% (the SIE lens for DR5) and thus are consistent with 
their result. 

Next, we calculate the p-q contours using the SIS lens 
and the SIE lens for SDSS DR1 and DR5 data. We limit 
our investigation to p < 5 and q < 3 because the value of 
4>(a; z) become unphysical beyond this range. The reason 
is as follows. We consider redshift evolution of the number 
density of lensing galaxies by multiplying current velocity 
function, ©, by the velocity function ratio, 4>(a; z)/<fi(a; z = 
0) (see (2.4.2)) just like Mitchell et al. (2005). If we try to 



calculate the ratio in the case of p > 5, the resulting ratios 
become negative. Moreover the parameter a in the velocity 
function becomes negative in the case of q > 3. 

Our results of the likelihood functions of the parame- 
ters, (p, q) , are in Fig. 2. In Fig. 2, we obtain maximum val- 
ues in L at (p = 5.0, q = 0.1) for the SIS lens and (p = 4.5, 
q = 0.1) for the SIE lens which are largely different from 
the estimation from the simulation (jp — 0.3, q — 0.75 the 
white-dots in Fig. 2). 

The maximum value, p and q are at the very corner 
of the parameter space so we expect that our best fit val- 
ues of p and q are the maximum of the likelihood within 
our physically meaningful limits and it is probably not a 
true maximum of the likelihood. But we think it still can be 
said that the values for p and q used by Mitchell et al. are 
inconsistent with the presently available observational data. 

Furthermore we compare the power-law evolution 
model and the model given by Mitchell et al. (p = 0.3, 
q = 0.75). We use our best fit parameter value of the power- 
law evolution model, v v , in the case of SIE lens, (y v — —0.15 
for DR1 and v v — —0.25 for DR5) and calculate the ratio 
of the velocity function at z — 0.3 and z — 1.0 to that at 
z — 0.0. The error of the parameter v n is large, so we use 
only estimated v v and set v n — 0. Similarly, we use the 
Mitchell et al.'s evolution model(p = 0.3, q — 0.75)and cal- 
culate the ratio of the velocity function. This result is shown 
in Fig. 3. As a result, we find that the evolution model by the 
Mitchell et al. has possibility to underestimate the number 
density of galaxy whose velocity dispersion is large. 



5 CONCLUSIONS AND DISCUSSION 

We employed the 'lens-redshift test' proposed by Kochanek 
(1992) and constrain the parameters of two evolution mod- 
els, the power-low evolution model and the evolution model 
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Figure 1. The likelihood functions of the number density evolution index u n and the velocity dispersion evolution index u v . The functions 
are normalized the corresponding maximum (Lmax)- The cross represent the maximum likelihood values. Top left panel: We use the SIS 
lens and the DR1 data as the parameters of the velocity function. The maximum value of L is u n = —2.0, v v = 0.15. Top right panel: We 
use the SIS lens and the DR5 data as the parameters of the velocity function. The maximum value of L is u n = —1.0, v v = 0.0. Bottom 
left panel: We use the SIE lens and the DR1 data as the parameters of the velocity function. The maximum value of L is v n = 1.5, 
h> v = —0.15. Bottom right panel: We use the SIE lens and the DR5 data as the parameters of the velocity function. The maximum value 
of L is v n = 2.0, v v = —0.25. 



proposed by Mitchell et al. (2005). We use the SDSS lens 
sample that is well-defined. It turns out that one of the pa- 
rameters of the power-law evolution model, v v is consistent 
with that given by Ofek, Rix & Maoz (2003) and Chae & 
Mao (2003) and the other parameter v n is consistent with 
the result by Ofek, Rix & Maoz (2003). 

On the other hand, the most likely values for the pa- 
rameters used in the evolution model by Mitchell et al. 
(2005) are inconsistent with the predicted values from the 
presently available observational data. The inconsistency 
may be mainly caused by ignoring the total lensing rate 
in the present analysis. Although Chae & Mao (2003) have 
shown that the evolution model can be restricted more 
strictly by including the information of lensing rate, we have 
not included the effect because the total rate depends on the 
magnification bias which can not be calculated at the mo- 
ment with very good accuracy in a consistent way with our 
SDSS sample. Furthermore we find that there is a possibility 
to underestimate the number density of galaxy which have 
large velocity dispersion by using the parameters in Mitchell 
et al.'s evolution model. 

The SQLS is still ongoing and the number of the strong 
lensing events is steadily increasing. Thus if we are able to 



calculate accurate magnification bias based on SDSS sample, 
our method will be very useful to constrain the evolution 
models of the number density of the lensing galaxies and 
thus the CDM based structure formation scenario. 
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Figure 2. The likelihood functions of the parameters, (p, q) used in Mitchell's evolution model. The functions are normalized the 
corresponding maximum (L m ax)- Top left panel: We use the SIS lens and the DR1 data as the parameters of the velocity function. The 
maximum value of L is P = 5.0, Q = 0.1. Top right panel: We use the SIS lens and the DR5 data as the parameters of the velocity 
function. The maximum value of L is P = 5.0, Q = 0.1. Bottom left panel: We use the SIE lens and the DR1 data as the parameters 
of the velocity function. The maximum value of L is P = 4.5, Q = 0.1. Bottom right panel: We use the SIE lens and the DR5 data 
as the parameters of the velocity function. The maximum value of L is P = 4.5, Q = 0.1. The white points represent estimation from 
simulation. 
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Figure 3. Ratio of the velocity function at z = 0.3 and z = 1.0 to that at z = 0. The set, (p = 0.3, q = 0.75) is estimated by numerical 
simulated values. And DR1 and DR5 in the figure show the ratio of velocity function when we use the best fit values estimated by using 
SIE lens model, (y n = 0.0, u v = —0.15) and (y n = 0.0, u v = —0.25), respectively, (see figure 1 and section 4) 
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